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This  paper  is  the  first  in  a  series  to  analyse  the  ciccnracy  of  polynomial  in¬ 
terpolation  of  functions  and  its  dependence  on  the  locations  of  the  interpolation 
nodes.  It  surveys  known  results  for  polynomial  interpolation  -n  an  interval.  It 
also  introduces  the  concept  of  the  Minimal  Interpolation  Sets  which  are  the 
“optimal”  interpolation  sets.  New  results  concerning  the  properties  of  the  min¬ 
imal  sets  as  well  as  procedures  for  locating  the  minimal  sets  are  presented. 
The  table  for  the  minimal  sets  in  the  norm  is  given.  An  zwlaptive  scheme 
for  determining  the  interpolation  order  is  also  presented.  Examples  show  the 
efficacy  of  this  approach. 
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I.  Introduction 


Problem  of  polynomial  interpolation  of  functions  is  a  classical  topic  in  approxima¬ 
tion  theory.  It  is  well  known  that  the  positions  of  the  interpolation  nodes  influence 
drastically  the  quality  of  the  approximation  for  high  degree  polynomials.  Various 
theorems  addressing  this  issue  are  known  for  interpolation  in  an  one  dimensional  in¬ 
terval,  particularly  through  the  work  of  P.L.  Chebyshev.  In  contrast,  in  two  and  three 
dimensions  when  the  domain  of  interpolation  is  not  a  rectangle  or  a  cube,  little  is 
known  concerning  the  approximation  accuracy  of  the  interpolation. 

Accurate  polynomial  approximation  of  functions  is  critical  for  the  success  of  Finite 
Element  computation,  especially  for  the  p  and  h-p  versions  of  the  Finite  Element 
method.  Functions  as  well  as  shape  informations  can  thus  be  specified  from  their 
values  at  certain  discrete  points  which  can  be  selected  a  priori. 

The  present  paper  is  the  first  in  a  series  (see  also  [2])  to  address  the  question  of 
the  dependence  of  the  approximation  accuracy  on  the  interpolation  nodes.  It  gives 
a  survey  of  known  results  related  to  polynomial  interpolation  in  an  interval.  It  also 
introduces  the  concept  of  minimcil  interpolation  sets  and  presents  the  principle  behind 
the  computational  procedure  used  in  finding  these  minimail  interpolation  sets  for  an 
interval  as  well  as  for  a  triangle.  Finding  the  minimal  interpolation  sets  for  the 
triangle  is  the  focus  of  the  second  paper. 

In  section  II,  we  review  known  results  pertinent  to  our  study.  Section  III  intro¬ 
duces  the  concept  of  the  minimal  interpolation  sets  for  various  norms.  Interpolation 
using  the  minimal  sets  leads  to  the  minimal  interpolation  error.  The  L°°  minimal 
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sets  are  listed  in  Table  1  for  polynomials  up  to  order  20.  Fig.  4  gives  the  graph  which 
characterizes  the  error  bound  for  such  interpolation.  Section  W  discusses  interpola¬ 
tion  error  in  the  derivative  of  the  function.  In  Section  V",  we  show  various  numerical 
interpolation  examples.  Section  M  presents  a  procedure  for  the  adaptive  selection  of 
the  degree  of  the  interpolation  polynomial  for  a  given  error  tolerance. 

II.  Overview  of  Known  Results 

Let  I  =  [—1,1],  C{I)  be  the  space  of  real  continuous  functions  on  I,  Vn  be  the 
space  of  algebraic  polynomials  of  degree  less  or  equal  to  n. 

Polynomial  interpolation  of  a  function  is  to  find  an  approximate  polynomial  which 
takes  the  value  of  the  function  at  some  given  points  in  /.  Such  polynomial  always 
exists. 

Theorem  1  Let  /(•)  €  C(I),  T  =  (to, T,,...,r„),  where  -1  <  tq  <  rj  <  ...  <  r„  <  1, 
then  there  is  a  unique  polynomial  p„(-)  in  Vn  such  that  p„(ri)  =  /(r,),  i  =  0, 1, ...,  n. 
p„{  )  is  given  by  Lagrange’s  formula 
Pnit)  =  Crm  :=  ELofirMt),  where  L,it)  = 

Ct  is  a  linear  projection  operator  which  maps  a  real  continuous  function  to  its 
corresponding  polynomial  interpolajit.  n  is  called  the  order  of  the  interpolation  T.  Ti 
are  called  nodes  of  the  interpolation  T. 

Theorem  2  (Weierstrass)  [7]  Every  continuous  function  in  C{I)  can  be  uniformly 
approximated  by  polynomials  to  any  degree  of  accuracy. 

The  question  arises  as  what  is  the  best  polynomial  approximant  in  Vn  to  a  given 
function.  We  define  the  least  deviation  d„(/)  ;=  inf{suptei\q(t)  —  /(0l>9(0  €  Vn}. 
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This  least  deviation  is  always  attained. 

Theorem  3  [6]  Given  /(•)  €  C(I),  there  exists  )  €  Vn,  such  that  for  all  q„(  )  € 
'Pn,  11/  -  u^nlloo  <  11/  -  qnWoc.  where  II  •  Hoc  is  the  I*  norm  in  C(/),  i.e.,  H/H,,,  = 
5«Pfe/|/(0l- 

Let  uj{S,f{-))  be  the  modulus  of  continuity  of  f{t)  in  /,  i.e.,  u;(6, /{•))  ;= 
sup|,,_,ji«,,,,,,€/|/{<i)  -  /(t2)l.  then 

Theorem  4  (Jackson)  [7]  Let  /(•)  be  an  r  times  continuous  differentiable  function, 
f  be  its  r(th)  derivative,  then  d„(/)  <  /*''^(-)),  where  A  <  ns/^j2. 

Therefore,  the  least  deviation  d„(/)  decreases  with  n  at  least  algebraically  to  zero 
with  power  r  if  f  has  r  continuous  derivative.  Moreover,  Bernstein  showed  that 

Theorem  5  (Bernstein)  [7]  /(•)  is  analytic  on  /  if  and  only  if  lim„_oo(dn(/))'^’'  <  1. 
/(■)  is  entire  if  and  only  if  /im„-,oo(<^n(/))'^"  =  0. 

Since  Lj  is  a  linear  projection  operator  from  C{I)  to  Vn,  we  define  the  U  norm 
of  Ct,  1  <  p  <  oc,  as  |1£tI1p  :=  sup/#o./€C(r)^^f^i]^-  where  the  U'  norm  of  a  function 
/  €  /  is  ll/lip  =  (//,  |/(t)|'’df)‘-''. 

From  Theorem  1,  we  have 

||£t||oc  =  A(r)  :=  s«p,e/i:Mt)|. 

k=0 

The  L°°  norm  A(T)  is  also  called  the  Lebesgue  constant  of  Ct- 
The  “mean  norm”  of  Ct  is  defined  as 

II<£t)II  :=(/' 

fc=0 
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Note  that 


||£r||p  <  li:rlp  ■.=  (f\j:\L,(t)\Ydty'^’>. 

k=Q 

We  shall  call  l^Crlp  the  L**  pseudonorm  or  the  in'*  norm  of  the  interpolation  operator 
Ct- 

Since  Cj  is  a  projection  operator,  HjCtIIp  >  2*/'’. 

We  emphasize  that  the  norm  of  Ct  depends  only  on  the  distribution  of  the 
interpolation  nodes  r, . 

How  well  can  Ct  f  approximate  /?  We  have  the  following  Lebesgue  inequ2dity. 
Theorem  6  For  all  /(•)  €  C(/),  1  <  p  <  oo,  ||/  -  £r/||p  <  (2‘/p  +  ||£r||p)d„(/). 

The  proof  uses  the  method  of  “intermediate  approximation” .  Note  that  the  n-th 
order  interpolant  of  the  best  approximating  polynomial  u;„  in  Theorem  3  is  itself, 
hence 

11/  -  ^t/I1p  =  11/  -  +  CT(Wn  -  /)llp  <  11/  -  U.-„llp  +  lKT(l^'n  -  /)l|p  <  11/  - 

+  ||£r||p)  =  (2‘/P  +  11-C7-||pK(/). 

Therefore,  the  interpolation  error  of  a  function  is  always  less  than  a  numerical 
factor,  which  does  not  depend  on  the  interpolation  function  but  depends  sensitively  on 
the  distribution  of  the  interpolation  nodes,  multiplied  by  the  least  deviation.  Note  for 
a  given  interpolation  T,  the  bound  ||£r||pd„(/)  for  ||£r{ti;„— /)|lp  is  optimsd,  i.e.,  there 
exists  a  nonpolynomial  continuous  function  /such  that  ll£xllp(i„(/)  =  1|£7{«;„  — /)l|p, 
where  is  the  best  approximating  polynomial  in  Theorem  3.  In  this  sense,  we  say 
the  interpolation  error  bound  in  Theorem  6  is  optimal. 
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We  shall  study  the  efficacy  of  the  interpolating  polynomials  as  the  number  of 
interpolation  nodes  increases.  We  shall  consider  an  infinite  triangle  array  consisting 
of  T',  where  T'  =  (tq'*,  {T'}  defines  an  interpolation  scheme. 

Given  /,  What  is  the  convergence  property  of  Cj./l 

Theorem  7  [7]  For  any  {  T'  },  H-Ct* Iloo  >  2;r“^fog(n  +  l)  +  6„,  where  6„  is  a  bounded 
sequence. 

The  right  hand  side  of  this  inequality  is  related  to  the  norm  of  the  Fourier  convo¬ 
lution  operator  given  by  the  Dirichlet  kernel. 

Theorem  8  (Erdos)  [3]  There  exists  a  positive  constant  c  such  that  for  any  {  T'  }, 
ll^r«||oo  >  27r"‘/o5(n  -h  1)  -  c. 

It  is  possible  to  deduce  the  following  result  from  Theorem  7  or  Theorem  8, 

Theorem  9  (Faber)  [7]  No  matter  what  are  the  nodes  of  an  interpolation  scheme  ( 
7*'  },  i  >  1,  there  is  a  continuous  function  /(•)  such  that  the  sequence  of  interpolating 
polynomials  relative  to  the  given  scheme  diverges  in  norm. 

It  is  important  to  realize  that  such  functions  are  non-differentiable  in  /,  in  fact, 
they  are  not  even  piecewise  differentiable.  On  the  other  hand,  for  a  given  continuous 
function,  it  is  always  possible  to  choose  an  interpolation  scheme  (depending  on  the 
function)  such  that  the  interpolating  polynomials  converge  to  the  function  in  L°° 
norm. 

To  illustrate  the  importance  of  the  choice  of  the  interpolation  nodes,  consider  the 
uniform  interpolation  scheme  =  (—1  -f-  2i/n,  j  =  0, 1, ..., n},  then 
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Figure  1;  as  a  function  of  n. 

Theorem  10  [7]  For  a  uniform  system  of  nodes  A(r„"„,/)  >  Ceip(n/2),  where 
C  >  0. 

The  exponential  divergence  of  is  apparent  in  Fig.  1.  Theorem  10  im¬ 

plies  that  ||^r.\,J|p  diverges  at  least  exponentially  for  p  >  1.  This  follows  from  the 
following  lemma. 

Lemma  1  [5]:  Let  u  be  a  polynomial  of  degree  n,  then  for  1  <  p  <  q  <  oo, 
Hull,  <  (2n)Vp-*/.||u||^. 

Take  q  =  oo  in  Lemma  1,  it  follows  that  for  an  n-th  order  interpolation  operator 
ll^rllp  >  \\CT\U{2nflr 

Note  that  from  theorem  6,  /  —  Cj-f  converges  in  norm  for  the  uniform  interpo¬ 
lation  scheme  if  d„(/)  decreases  faster  than  However,  this  is  not  the  case 

for  functions  with  only  finite  number  of  derivatives.  Even  for  analytic  (or  infinitely 
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Figure  2:  The  function  R(t)  =  (1  -f-  and  C-ri  R.  R. 

differentiable)  functions,  this  condition  may  not  be  satisfied.  .A.s  Runge  showed  that 
for  R{t)  =  (1  -|-25t'^)“‘,  jR  diverges  from  R  in  intervals  0.728  <  |t|  <  1.  This 
is  illustrated  in  Fig.  2  for  ^R  and  £7-13  ^R.  We  see  that  as  interpolation  order 
n  increases,  the  interpolant  £t",  ^Rif)  diverges  farther  away  from  R(t)  in  the  speci¬ 
fied  intervals.  In  this  sense,  the  uniform  interpolation  scheme  is  generally  considered 
“not  good” .  Bernstein  gave  a  similar  result  for  the  piecewise  differentiable  function 
B(t}  =  |f|,  He  showed  that  diverges  from  B{t)  in  intervals  0  <  |t|  <  1. 

For  T  consists  of  the  zeroes  of  the  ( n  -f- 1  )-th  order  Chebyshev  polynomials  = 
{— i  =  0, ...,  n},  we  have 

Theorem  11  (Chebyshev)  [6]  <  2Tr~^log[n  -f  1)  -I-  47r“'(l  -|-  logA). 

In  view  of  Theorem  8,  w'e  see  that  the  Chebyshev  interpolation  scheme  {  } 

is  almost  optimal. 
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Note  that  other  than  in  I*  norm,  under  an  appropriate  choice  of  interpola¬ 
tion  scheme  {  T'  },  £7.)  may  converge  to  /  for  all  /  €  C(I).  In  particular,  if 
we  choose  T  to  be  the  quadrature  points  for  the  Gauss- Legendre  quadrature,  i.e.. 
Tigg  :=  (xi,  •••1  Jn+i)i  where  Xk's  are  the  zeroes  of  the  (n  -f  l)-th  order  Legendre 
polynomial,  P„+i(t),  then 

Theorem  12  [6]  If /(•)  €  C{I),  I\f-Crijh  <  2V2d„{f)-  \\f  -  <  4d„(/): 

Similarly,  we  can  also  choose  T  to  be  the  quadrature  points  for  the  Gauss-Lobatto 

quadrature,  t.e.,  Tl^  =  (-l,x', . x;.  1),  where  are  the  zeroes  of  the  derivative 

of  the  n-th  order  Legendre  polynomial,  P'nit)-  The  same  inequalities  hold  for  /  - 
£r«^^/.  We  note  that  for  similar  inequalities  hold. 

The  nodes  for  the  Fejer  interpolation  scheme  consists  of  the  zeioes  of  the  integral 
of  the  Legendre  polynomials,  i.e.,  where  x['s  are  the  zeroes 

of  the  integral  of  the  n-th  order  Legendre  polynomial,  /I,  Pr,(x)dx.  Note  that  x|  = 

-T-rUi  =  l- 

Interpolation  schemes  using  (  r;,,  { Tu ),  { }.  { }  will  be  called 

Legendre,  Lobatto,  Fejer  and  Chebyshev  schemes  respectively. 

Maday  has  recently  shown  that  if  /(•)  €  //'"(/),  t.e.,  if  edl  the  derivatives  of  /  up 
to  order  m  are  square  integrable,  then  converges  to  /  in  if  m  >  1: 

Theorem  13  (Maday)  [1]  If  /  €  //"*(/),  then  there  is  a  positive  constant  c,  such 
that  11/  -  £7»^^/||,^i(/,  <  cTi‘“'"||/||ff-.(/). 

The  i/"*  norm  for  /  is  ||/||h>"  =  (fl,  where  d'f(t)  denotes  the 

i-th  derivative  of  /. 
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Erdos  has  shown  the  following  results  concerning  the  I*  pseudonorm  (|L|'  norm) 
and  the  mean  norm  for  the  interpolation  operator 

Theorem  14  (Erdos)  [3]  There  exists  a  positive  consiam  C  such  that  for  any  n-th 
order  interpolation,  |jCt|i  >  Clog(n  +  1). 

Theorem  15  (Erdos)  [3]  To  every  positive  e,  there  exists  a  positive  integer  .V  such 
that  for  every  ?i-th  order  interpolation  operator  T  satisfying  n  >  .V,  IK^r)!!'  >2  —  6 

By  definition  of  the  Gauss- Legendre  quadrature,  we  have  for  the  Legendre  inter¬ 
polation  scheme,  =  2. 

Theorem  16  (Fejer)  [4]  The  only  n-th  order  interpclation  set  for  which  |L*(0I^  < 
l,t  €  /  is  the  n-th  order  Fejer  set 

Note  that  and  include  the  end  points  {-1.1}  of  I.  Henc"',  if  a  function 

is  piecewise  approximated  Ly  polynomial  interpolation  in  contiguous  intervals  using 
the  Lobatto  sets  or  the  Fejer  sets,  its  continuity  at  the  joint  points  is  guaranteed. 
This  property  is  essential  for  the  continuity  of  the  piecewise  interpolated  polynomial 
function  in  an  interval. 

III.  Minimal  Interpolation  Sets 

In  the  following,  we  consider  interpolation  schemes  whose  nodes  include  the  end 
points  {-1,1}.  From  Theorem  6,  for  a  given  interpolation  T,  the  optimal  upper  bound 
for  the  interpolation  error  of  a  function  /  in  L’^  norm  is  (1  -f-  |j£7-l|oc)<^n(/)-  We 
can  minimize  this  optimal  error  upper  bound  by  minimizing  the  Lebesgue  constant 
of  the  interpolation  operator  K^tIIoo  through  the  redistribution  of  the  interpolation 
nodes  (the  least  deviation  d„{f)  by  definition  is  already  the  smallest  possible  error). 
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Similarly,  the  interpolation  error  of  a  function  /  in  norm  is  less  than  (2*'^  + 
(|£Tjlp)d„(/V  However,  11£t|1p  is  hard  to  work  with.  Since  |l£r!l.r  <  |£r|p.  (2'''"  + 
\^T\p)dn{})  also  gives  an  upper  bound  for  the  interpolation  error  in  norm,  we  can 
instead  minimize  |£tIp  which  has  a  closed  form.  We  shall  see  that  the  corresponding 
minimal  sets  are  also  good  interpolation  sets.  We  shall  say  that  an  interpolation  sei  is 
minimal  in  \L\^  or  Z,“  if  it  minimizes  |£t|p  or  ||£r||(»-  We  shall  say  an  interpolation 
set  is  minimal  in  the  mean  if  it  minimizes  |K£r)||-  order  minimal  sets  are  denoted 
as  Tjtip  and  respectively.  Interpolation  schemes  using  {  }.  {  } 

and  (  }  are  called  minimal  interpolation  schemes  in  \L\^,  I*  and  in  the  mean. 

We  shall  see  that  minimal  interpolation  sets  have  quite  similar  properties. 

i.  Minimal  interpolation  sets  in 

For  a  given  system  of  interpolation  nodes  T  =  (tq,  Tj,  ...,  r^),  where  tq  =  —l.Tn  = 
1,  let  A,(r)  =  then  A(r)  =  77iax,=, . „A.(j  ;. 

Conjecture  1  T  minimizes  ||£tI1oo  if  and  only  if  A|(r)  =  \2(T)  =  ...  =  A„(T),  these 
n  —  1  conditions  uniquely  determine  the  n  —  1  interpolation  nodes  Tj,  ..  ,  rn_i. 

Note  that  the  condition  that  the  minimal  sets  contain  the  end  points  {-1,1}  of  the 
interval  I  is  essential.  The  minimal  sets  are  not  unique  if  {-1,1 }  is  not  included  in  the 
minimal  sets  As  can  be  verified  in  the  case  of  order  2  interpolation,  any  interpolation 
set  7^  =  (-y,0,y)  where  2>/2/3  <  y  <  1  is  also  minimal. 

The  uniqueness  implies  that  the  minimeil  sets  are  symmetrical,  i.e.  they  are  in¬ 
variant  under  the  reflection  x  -*  —x. 
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Figure  3;  The  coordinates  for  T^l. 

We  have  numerically  found  the  sets  satisfying  the  equalities  in  Conjecture  1  for 
interpolation  order  up  to  n  =  50.  They  appear  to  be  minimal  in  .  We  shall  denote 
n-th  order  such  set  by  Conjecture  1  states  that 

In  the  appendix,  we  list  the  coordinates  for  T/**  from  n  =  1  to  n  =  20.  The  error 
in  the  coordinate  is  less  than  10~'.  Fig.  3  shows  the  coordinates  of  They  are 
almost  uniformly  distributed  in  the  circle  coordinate:  0  =  — cos~'(i), 

Asymptotically  AIT/**,)  scales  with  login  +  1).  Fig.  4  shows  the  Lebesgue  con¬ 
stant  (the  L°°  norm)  for  T/*®.  Our  numerical  result  gives  asymptotically  A{Tt«,) 
27r"'/og(n  +  l)-t-0.522  ~  0.523.  This  is  consistent  w'ith  Theorem  8.  We  note  that  even 
at  order  50  (which  is  considerably  higher  than  the  interpolation  order  in  any  practical 
situation),  A{T,-»)  is  only  approximately  3,  hence,  from  Theorem  6,  the  interpolation 
error  would  be  at  most  4  times  the  least  deviation  at  order  50.  Note  the  Lebesgue 
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Figure  4:  A(Tt*oe)  as  a  function  of  n. 

constant  of  the  interpolation  operator  depends  very  sensitively  on  the  distribution 
of  the  interpolation  nodes.  For  the  uniform  distribution,  from  Fig.  1,  the  Lebesgue 
const  emt  at  order  20  already  exceeds  10'*. 

The  |Z,|'’  norm  of  also  appears  to  scale  with  log(n  +  1).  However,  the  square 
of  the  mean  norm  saturates  from  below  at  2. 

ii.  Minimal  interpolation  sets  in  |L|',  |I|^  and  in  the  mean 

The  minimal  interpolation  sets  in  |£|',  \L^  and  in  the  mean  are  found  by  using 
some  standard  minimization  procedure  for  the  corresponding  norms. 

We  shall  illustrate  in  this  section  that  T|2|i,  T^l^2,  have  similar  scaling  be¬ 
haviors  as  T^=o.  Their  coordinate  positions  are  relatively  close,  and  their  \L\^  norms, 
which  determine  the  upper  bounds  for  the  interpolation  error,  are  almost  the  same. 

It  was  conjectured  [3]  that  the  n-th  order  minimal  set  for  the  mean  norm  is 
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Figure  5:  the  differences  in  the  Lebesgue  constant  (the  L**  norm)  between  Tw  and 
(diamonds),  7]^|a  and  (triangles),  and  (asterisks),  and  T/** 
(squares)  respectively. 

the  n-th  order  Fejer  set.  The  mean  minimal  sets  we  have  found  indicate  that  this 
conjecture  is  incorrect. 

Fig.  5  shows  the  differences  in  the  Lebesgue  constant  (L°°  norm)  between  T|]^|, 
and  T[^oo,  T”^2  and  T);"*,  and  and  respectively.  These  differences 

appear  to  asymptote  to  constants.  In  the  case  of  and  Tj^ja,  it  appears  that 
they  are  very  close  to  T^"®  as  n  increases  (the  asymptotic  differences  in  the  Lebesgue 
constants  are  close  to  zero).  The  differences  in  the  Lebesgue  constants  between  all 
these  sets  are  small  so  that  these  sets  are  almost  the  same  in  practical  interpolation 
situations.  E.g.,  if  we  use  instead  of  7^"®  the  difference  in  the  Lebesgue  constant 
(L°°  norm)  is  only  of  order  0.07.  According  to  theorem  6,  this  only  increases  the 
interpolation  error  for  approximately  0.07  times  the  least  deviation,  which  is  usually 
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Figure  6:  The  mean  norm  for  the  mean  minimal  set  ||(£r» .  )||^  as  a  function  of  n. 
negligible  in  practice. 

The  norm  as  well  as  the  \L\^  norm  for  all  minimal  sets  scale  with  log{n  +  1). 
The  square  of  the  mean  norm,  ||(£t)||^.  asymptotes  to  2  from  below  for  all  minimal 
sets.  Fig.  6  shows  the  dependence  of  the  mean  norm  for  the  mean  minimal  set 
on  the  interpolation  order  n. 

iii.  The  Lobatto  sets 

For  the  Lobatto  sets,  using  the  Gauss-Lobatto  quadrature  formula,  one  obtains 
that  the  square  of  the  mean  norm  equals  to  2  —  C„,  where  C„  >  0  and 

/im„_ooC’n  =  0.  From  Theorem  15,  one  sees  that  the  Lobatto  sets  are  almost  optimal 
in  the  mean  norm. 

The  Lobatto  sets  are  usually  good  enough  for  interpolation  purpose  in  practice. 
(As  is  any  set  consisting  of  the  zeroes  of  any  other  orthogonal  polynomial  in  I.) 
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However,  albeit  there  are  generalization  of  orthogonal  polynomials  in  some  standard 
connected  domains  in  several  dimensions,  their  relations  to  quadrature  and  interpo¬ 
lation  are  unclear.  The  minimal  sets,  on  the  other  hand,  have  direct  generalization  in 
these  situations  and  in  view  of  Theorem  6,  they  are  the  “optimal”  interpolation  sets 
[2]. 

IV.  Approximation  of  Derivatives  of  Functions 

Frequently,  we  not  only  need  to  approximate  the  function,  but  also  its  derivative. 

Consider  the  following  problem.  Suppose  we  have  only  access  to  the  function  val¬ 
ues,  but  not  its  derivative  values  at  certain  points,  we  want  to  extrapolate  information 
about  its  derivative.  Theorem  1  again  provides  a  way  to  estimate  the  derivative  of 
the  function  by  evaluating  the  derivative  of  the  interpolating  polynomial; 

p'„(0  =  £V/(0 

Here  prime  (^)  denotes  taking  derivative  with  respect  to  t.  We  call  the  derivative 
interpolation  operator.  We  note  that  it  b  not  a  projection  operator. 

Let  g„  be  an  n-th  order  polynomial,  using  the  method  of  intermediate  approxima- 

lion,  II/'  -  Cr/t  <  -  S'.IU  +  IICtIIpII/  - 

||£i-IU  =  V(r)  := 

114-11.  <  141,  :=  ( /' (t  liifOIW"- 

fc=0 

Again  we  shall  call  \C''p\p  the  V  pseudonorm  or  the  \[/\  norm  of  the  derivative 
operator  and  A'(7’)  the  Lebesgue  constant  for  the  derivative  operator. 
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Figure  7;  The  Lebesgue  constants  for  the  derivative  interpolation  operator  for 
(triangles),  (squares),  (asterisks)  and  T^,,j(diainonds). 

The  approach  we  have  used  in  the  previous  section  can  be  used  to  minimize  \C'y\p. 
We  shall  minimize  |£V1p  mean  norm  for  the  derivative  interpolation  operator 

ll(^r)ll  •=  (f-i  S2=o  .  We  shall  denote  the  corresponding  minimal  sets 

by  and  respectively. 

The  error  for  the  interpolation  of  the  derivative  is  usually  largest  at  the  end 
points.  In  quite  contrast  :o  i  hf  Lebesgue  constant  for  the  interpolation  operator,  the 
Lebesgue  constant  for  the  derivative  interpolation  operator  A'(T)  appears  to  increase 
quadratically  with  n  for  both  and  T^q,.  Fig.  7  shows  this  scaling  of  the  Lebesgue 
constant  for  the  derivative  interpolation  operator  for  Tj^,|j,  T^.),  T£oo  and 

However,  the  |X.|'’  and  the  mean  norm  for  the  derivative  interpolation  operator 
appear  to  increase  linearly  with  the  interpolation  order  n.  As  is  illustrated  in  Fig.  8 
for  the  lip  norm.  This  is  a  result  similar  to  Theorem  13. 
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Figure  8;  The  \L\^  norms  for  the  derivative  interpolation  operator  for  (triangles), 
squares),  (asterisks)  and  rto4(diamonds). 

We  note  that  the  values  for  the  Lebesgue  consteints  for  the  derivative  interpolation 
operator  for  different  minimal  sets  are  quite  close.  This  indicates  that  the  minimal 
sets  for  the  derivative  interpolation  operator  are  quite  close  to  the  minimal  sets  for 
the  interpolation  operator. 

Suppose  we  are  not  only  able  to  access  the  function  value,  but  also  its  derivative 
at  certain  n+1  points.  Then  we  can  use  Hermite  interpolation  formula  to  obtain  a 
2n  +  1  order  polynomial  to  approximate  the  function; 

P2n+i(f)  =  :=  T.l=o{f(Tk)Hk{t)  +  f'(Tk)Hk{t)),  where  Hk{t)  =  (1  - 

2L'k(Tk){t  -  Tk))Ll{t),  Hkit)  =  {t-  Tk)Ll{t). 

The  minimal  sets  for  the  Hermite  interpolation  operator  Hr  can  be  considered 
similarly. 
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V.  Examples 


Let  us  consider  the  interpolation  of  the  following  rational  function  =  7:^^- 
In  Fig.  9,  we  show  the  interpolation  error  Ri{t)  —  CrRiit)  for  minimal  sets 

^(t)  Lobatto  set  for  order  n  =  9  interpolation.  Interpolation 

using  minimal  sets  does  appear  to  have  smaller  interpolation  error  than  any  known 
interpolation  schemes.  The  magnitudes  of  the  interpolation  error  for  dilferent  minimal 
sets  are  of  the  same  order.  Indeed,  T/l*,  appears  to  have  the  smallest  interpolation 
error  and  has  the  largest  error  while  the  interpolation  error  for  seems  to  lie 
in  between.  The  interpolation  errors  for  T|^|j  and  appear  to  be  very  close  to  that 
of  T^cc.  This  is  consistent  with  the  fact  that  the  norms  for  T^i^2  and  7y[|i  are  very 
close  to  the  norm  of  T/L,.  The  ordering  of  the  norms  for  minimjd  sets  appear  to 
be  indicative  of  the  ordering  of  their  interpolation  errors. 

On  the  other  hand,  if  we  interpolate  the  same  function  to  extrapolate  the  deriva¬ 
tive,  the  error  is  larger.  We  show  in  Fig.  10  the  derivative  error  —  CrRiit) 
for  minimal  sets  T(V)'  ^06  fo**  interpolation  order  n  =  9.  The  error 

in  derivative  is  largest  at  the  end  points  f  =  ±1.  T^n  appears  to  have  the  smallest 
derivative  error,  while  T^«,  has  the  largest  error.  They  all  have  the  same  order  in 
interpolation  error. 

For  minimal  sets,  both  the  interpolation  error  in  norm  or  in  H*  norm  for 
/?/3  decreases  exponentiaJly  (see  Fig.  11)  ,  in  accordance  with  Theorem  12  and  the 
fact  that  the  function  R$  is  analytic  hence  d„{Rfi)  decreases  exponentially.  The 
interpolation  error  in  norm  is  dominated  by  the  error  in  its  derivative,  since  as 
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Figure  9:  The  interpolation  error  Ri(t)  —  CtRiH)  for  minimal  sets  T|'}_|2  (dot),  T|||i 
(dash  dot),  T£*®(solid),  (dash)and  the  Lobatto  set  (long  dash)  for  order 


Figure  10:  The  derivative  error  R'i{t)  —  C'rRiit)  for  minimal  sets  (dash), 
(dot),  T^oo  (solid)  and  (dash  dot)  for  interpolation  order  n  =  9. 
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Figure  11:  The  interpolation  error  in  norm  (11a)  and  norm  (11b)  for  where 
3  =  1  and  3  =  3. 

we  have  seen,  the  error  in  the  derivative  of  the  function  from  interpolation  is  much 
larger  than  the  error  in  the  function  itself.  All  minimal  sets  seems  to  have  almost  the 
same  order  of  interpolation  error  in  and  norms. 

VI.  Adaptive  Determination  of  Approximation  Order 

Suppose  we  want  to  approximate  a  function  in  an  interval.  It  is  neither  appropri¬ 
ate  nor  practical  to  use  polynomial  interpolation  for  the  function  in  the  entire  interval. 
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Figure  12:  The  polynomial  interpolation  of  the  function  /  =  in  4  subintervals 

/,  =  [0, 2],  h  =  [2, 4],  h  =  [4, 6],  h  =  [6, 8) 

Approximation  of  functions  using  polynomials  of  arbitrarily  large  order  may  be  pro¬ 
hibitively  expensive.  This  problem  could  be  solved  by  dividing  the  whole  interval 
into  subintervals  and  by  successively  approximating  the  function  in  each  individual 
subinterval.  Note  that  if  we  use  the  minimal  sets  discussed  in  section  III  which  con¬ 
tains  the  end  points  of  the  interval,  the  continuity  of  the  function  at  junction  points 
is  automatically  guaranteed. 

However,  if  we  use  the  same  polynomial  order  for  interpolation  in  all  the  subinter¬ 
vals,  the  approximation  error  may  be  of  quite  different  order  of  magnitude  in  different 
subintervals.  For  example,  in  Fig.  12.  we  show  the  polynomial  interpolation  of  the 
function  /  =  (/-|0)a+i  4  subintervals  /j  =  [0,2),/2  =  [2, 4], /a  =  [4,6],/4  =  [6,8] 

using  the  same  order  9  polynomial  interpolation.  We  see  that  the  function  is  well 
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Figure  13;  The  error  /  —  Crf  in  different  subintervals  for  order  9  interpolation. 

approximated  in  each  subinterval,  however,  as  is  shown  in  Fig.  13,  the  error  f  —  Crf 
varies  with  quite  different  order  of  magnitude  in  different  subintervals.  In  fact,  com¬ 
pared  to  the  error  in  U,  errors  in  other  subintervals  are  negligible.  Frequently,  it 
is  inefficient  to  have  such  a  situation,  since  both  the  access  of  the  function  values 
at  interpolation  nodes  and  the  interpolation  procedure  cost  more  for  high  order  in¬ 
terpolation  than  for  low  order  interpolation.  The  interpolation  order  in  subintervals 
h,h^h  can  in  fact  be  reduced  to  achieve  the  same  order  of  accuracy  (10~®  in  all 
subintervals).  As  is  illustrated  in  Fig.  14,  to  achieve  the  same  order  of  accuracy  for 
the  approximation,  we  only  need  order  5  interpolation  in  Ii ,  order  6  in  I2,  and  order 
7  in  I3  . 

In  the  following,  we  give  a  simple  adaptive  procedure  for  the  determination  of  the 
approximation  order.  Let  {  T'  =  =  0.  ...,n, ...;  4'^  =  — 1,t;^'^  = 

1,  =  0  }  be  the  interpolation  scheme  in  /  =  [— 1, 1].  Since  low  order  interpolation 
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Figure  14;  The  error  /  -  £7/  in  different  subintervals  for  interpolati-  n  using  adap¬ 
tively  determined  order. 


informations  are  relatively  cheap  to  obtain,  the  interpolation  error  can  be  estimated 
by  using  low  order  quantities.  At  order  n,  we  have  at  our  disposal  queintities  { 

=  0.  =  1 . n  }.  Define  the  estimated  interpolation  error  at  order  n 

l>y  =  maXj=o . .:.=o . /)( '^ )|.  Fot  a  prescribed  error  toler¬ 

ance  p(/),  calculate  the  estimated  interpolation  error  e„(/).  If  e„(/)  <  p{f).  stop; 
otherwise,  continue  to  go  to  higher  order  polynomial  interpolation  until  e„(/)  <  p(/). 


Fig.  14  is  generated  using  this  adaptive  procedure.  We  see  that  it  consistently 
gives  the  appropriate  interpolation  order  for  different  subinterv^lls. 

Numerical  computation  indicates  that  for  all  the  minimal  sets  we  have  consid¬ 
ered  as  w'ell  as  the  Lobatto  interpolation  schemes,  the  above  error  estimator  satisfies 


jEdn. 


1,  j.e.,  e„(/)  is  an  asymptotic  error  estimator  for  these  sets.  This  is 


ll/-£r»/lU 
because  all  these  sets  are  dense  in  I 
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Appendix 


The  coordinates  of  Minimal  Sets.  Minimal  sets  are  symmetrical  under  the  re¬ 
flection  X  —X.  Only  nonnegative  coordinates  are  given, 
n  =  1 

1.000000000000000 
n  =  2 

0.000000000000000  1 .000000000000000 
n  =  3 

0.4177912974318541  1.000000000000000 
n  =  4 

0.0000000000000000  0.6209113447507076 
1.000000000000000 
n  =  5 

0.2689070754306328  0.7341271205122524 
1.000000000000000 
n  =  6 

0.0000000000000000  0.4461215875027044 
0 . 8034402333430620  1  .OOOOOOOOOOOOOOO 
n  =  7 

0.1992876906268561  0.5674310617087234 
0.8488719620155624  1.000000000000000 
n  =  8 
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0.0000000000000000  0.3477880851677845 
0.6535334678686619  0.8802308501922200 
1.000000000000000 
n  =  9 

0.1585652002417824  0.4601504943340195 
0.7166138903070984  0.9027709837119849 
1.000000000000000 
n  =  10 

0.0000000000000000  0.2848881583261380 
0.5466676817797121  0.7640984500324875 
0.9195087505000881  1.000000000000000 
n  =  11 

0.1317517498523124  0.3862689500044784 
0.6144355686309110  0.8006822718709737 
0.9322747835437754  1 .000000000000000 
n  =  12 

0.0000000000000000  0.2412238223627454 
0.4684174854953188  0.6683666122414075 
0.8294354770322080  0.9422316270662717 
1.000000000000000 
n  =  13 

0.1127326359341456  0.3325422077453501 
0.53566551854944670.7119103291764621 
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0.8524275965461744  0.9501460629121580 
1.000000000000000 
n  =  14 

0.0000000000000000  0.2091512677285655 
0.4091565197775405  0.5912705391349697 
0.7475281136999568  0.3710916050399280 
0.9565402629254455  1 .000000000000000 
n  =  15 

0.0985297784103952  0.2918018772197364 
0.4738547252988244  0.6376896893581207 
0.7770061975121303  0.8864437448634451 
0.9617797393250818  1.000000000000000 
n  =  16 

0.0000000000000000  0.1845994015541723 
0 .362909644755691 4  0.528857282581 0849 
0.6767882746653493  0.8016617880341823 
0.8992200395501156  0.9661264747628877 
1.000000000000000 
n  =  17 

0.0875146300475834  0.2598845240957428 
0.4243549140109628  0.5759276777085162 
0.7099951825658469  0.8224813031663734 
0.9099637721447408  0.9697722157216912 
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1.000000000000000 
n  =  18 

0.0000000000000000  0. 1 652022750 114381 
0.3258963792565789  0.4776989264800610 
0.6164674647713962  0.7384152647561639 
0.8402138563369347  0.9190827135791494 
0.9728598817268702  1.000000000000000 
n  =  19 

0.0787199876854192  0.2342217468993406 
0.3839541966628281  0.5242305008694372 
0.6515953463715364  0.76291 13934252203 
0.8554359783038173  0.9268876580527369 
0 .975497771 24 1 1 774  1 .000000000000000 
n  =  20 

0.0000000000000000  0.1 4949 1 70707381 74 
0. 2956424244330606  0.435 1 88 1 757862846 
0.5650105307170965  0.6822087048389735 
0.7841637256141400  0.8685966695651883 
0.9336186965721600  0.9777691371910347 
1.000000000000000 
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